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The prevalence of Plasmodium falciparum malaria in Zanzibar has reached historic lows. Improving control 
requires quantifying malaria importation rates, identifying high-risk travelers, and assessing onwards 
transmission. 

Estimates of Zanzibar's importation rate were calculated through two independent methodologies. First, 
mobile phone usage data and ferry traffic between Zanzibar and mainland Tanzania were re-analyzed using 
a model of heterogeneous travel risk. Second, a dynamic mathematical model of importation and 
transmission rates was used. 

Zanzibar residents traveling to malaria endemic regions were estimated to contribute 1-15 times more 
imported cases than infected visitors. The malaria importation rate was estimated to be 1.6 incoming 
infections per 1,000 inhabitants per year. Local transmission was estimated too low to sustain transmission 
in most places. 

Malaria infections in Zanzibar largely result from imported malaria and subsequent transmission. 
Plasmodium falciparum malaria elimination appears feasible by implementing control measures based 
on detecting imported malaria cases and controlling onward transmission. 



The expansion of road, rail and air transport networks and their increasing affordability contribute to the 
spread of vector-borne diseases through human and vector movement^' ^. The relevant concept for malaria 
was named "vulnerability" by the World Health Organization (WHO) and defined qualitatively as ''the 
frequent influx of infected individuals or groups and/ or infective anophelines"^; quantitatively, the rate of malaria 
importation includes all parasites that cross the border in humans and vectors. The impact of vulnerability 
depends on an area's "receptivity" to malaria which reflects the conditions of transmission ''through the abundant 
presence of vector anophelines and the existence of other ecological and climatic factors'^ . Receptivity is defined 
quantitatively as the effective reproduction number R^, which describes the expected number of secondary human 
infections originating from a single, untreated infected human taking into account vector control measures^' ^. As 
part of renewed effort towards malaria elimination addressed in the Global Malaria Action Plan (GMAP), the 
malaria Eradication Research Agenda (MalERA) highlighted the importance of modeling strategies in defining 
vulnerability to plan for sustainable elimination and prevent reintroduction^' ^, especially in areas where malaria 
prevalence has been reduced but receptivity is still favorable to transmission^. 

In Zanzibar where the main malaria pathogen is Plasmodium falciparum, the introduction of artemisinin- 
based combination therapy in 2003 and the escalating ownership of insecticide-treated nets^' ^° resulted in an 
important decrease in malaria prevalence (Figure Therefore, the Zanzibar Minister of Health commissioned 
a report to assess the feasibility of malaria elimination^^' ^^ which highlighted the risk that cases may be imported 
from surrounding malaria-endemic regions of mainland Tanzania (Figure 1). Recent approaches to quantifying 
human mobility patterns point to novel insights through the analysis of mobile phone records and initial 
estimates of malaria importation based on such records estimated between 1 and 12 infections imported each 
year per 1,000 Zanzibar! citizens^^. Uncertainty in the estimated importation rate spanned a critical value; at the 
lowest estimates of malaria importation rates, malaria elimination was cost- saving and sustainable, but at the 
upper estimates, controlling malaria was optimal. 
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Figure 1 | Map of Zanzibar, (a) Zantel coverage regions in mainland Tanzania^^, (b) Zanzibar and its prevalence rate (%) on the two main islands: 
Unguja and Pemba 



Strategic planning for the future of malaria control therefore 
demands refining the estimate of malaria importation rates to nar- 
row the uncertainty, which could sway the decision of when and 
whether to eliminate malaria compared with a strategy of control 
in Zanzibar or elsewhere in Africa^^. Here, the concept of vulnerabi- 
lity has been rigorously defined and analyzed to identify high-risk 
travelers and the most important paths of malaria importation. In 
addition, mathematical models assuming heterogeneous exposure to 
malaria, recent malaria surveys, and mobile phone usage records and 
ferry traffic data were used to re-estimate malaria importation rates, 
cross-validate those estimates, and narrow the uncertainty. Finally, 
the interplay between malaria importation and subsequent local 
transmission was evaluated through the relationship between vul- 
nerability, receptivity and local malaria dynamics. 

Results 

The islands of Zanzibar have strong transport connections to main- 
land Tanzania where transmission levels are higher. While daily 
flights primarily bring tourists from non- endemic regions, ferry tra- 
vel from Dar es Salaam and informal movement on small fishing 
and trading vessels likely represent the highest risk pathways for 
imported infections. 

Contributions of residents, visitors and immigrants to imported 
malaria. A recently published mathematical model identified four 
ways that malaria parasites can be imported into a region^^: human 
residents of the region can become infected while traveling abroad 
and then return home for the duration of the infection; human 
visitors can spend time in Zanzibar after being infected elsewhere; 
human immigrants infected elsewhere can move permanently into 
the region; and infected mosquitoes can travel into the region 



(Figure 2). The frequency of importation through each one of 
these routes differs, as does the significance of each importation 
event for potential onwards transmission; some infections are 
unlikely to contribute to local transmission depending on the 
receptivity of the place where infections were imported and the 
duration of stay. Each untreated malaria infection imported in a 
resident or an immigrant would be counted the same; both are 
likely to remain in the country as long as they are infectious, 
providing a long window during which they may infect mosquitoes 
locally, even though the number of immigrants is likely very smalP^. 
Each infected human visiting the region, however, tends to have 
short visits compared with the long infectious period of malaria. 
To produce a standardized importation rate, each one of these 
visitors must be discounted by the length of stay relative to the 
duration of an untreated infection^^. Since Zanzibar is an 
archipelago separated from the mainland by approximately 40 km, 
malaria infections imported by mosquitoes are assumed to be rare, 
and they have not been considered further here^^. 

Analysis of length-of-stay from phone records in Tanzania sug- 
gests that the median visitor trip length to Zanzibar was approxi- 
mately 13 days, while the duration of untreated malaria infections is 
on average 200 days^^ If there were roughly as many infected return- 
ing residents as there were visitors to a region, then on average, each 
resident who returned infected and who was not properly diagnosed 
and treated would count for 15 infected visitors. The relative import- 
ance of transmission routes can be checked by doing a sensitivity 
analysis on the number of exogenous infections brought by residents 
and visitors and the length of stay of visitors relative to the duration 
of infection (Figure 3). For example, with similar numbers of 
exogenous infections and visitors staying 10% of the infection dura- 
tion, residents contribute to 90-95% of imported malaria. These 
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Figure 2 | Qualitative assessment of malaria importation to Zanzibar and its influence on local malaria dynamics. The pathogen can be introduced into 
the island according to four ways: 1) Zanzibar residents (R) can become infected while traveling abroad and bring it home; 2) infected visitors (V) to 
Zanzibar can bring the pathogen from mainland Tanzania; 3) infected immigrants (M) travel and settle permanently in Zanzibar; 4) infected vectors 
travel to Zanzibar either flying or passively carried out by wind or riding along human means of transportations. Once introduced, the probability for 
further transmission is based on the length and the receptivity of the area where the pathogen is imported, defining the net importation rate. 



estimations assume both residents and visitors stay in the same 
places in Zanzibar. 

Estimating the vulnerability of Zanzibar. Simply counting the 
number of travelers to Zanzibar may not provide accurate 
estimates of malaria importation, so estimation of malaria risk was 
done in two ways for cross-validation. First, estimates of importation 
were made from mobfle phone usage records and ferry traffic data. 
The proportion of residents infected whfle travelling abroad can be 
estimated from the sum of their time spent at various levels of risk. 
An important caveat is that infection risk whfle traveling is highly 
heterogeneous^^' ^\ Second, a mathematical model was developed 
that could explicitly take into account the dynamics of trans- 
mission and malaria importation; this model was apphed to cross- 
sectional parasite surveys from Zanzibar in 2007^^ to estimate the 
importation rate and R^. 

The first method, based on mobfle phone usage data, ferry traffic 
and heterogeneous risk of infection gives an estimate of between 0.7 
and 6 infections imported by residents per 1,000 Zanzibaris, based on 
the risk of malaria for the entire regions of origin, or their principal 
cities, respectively (Table 1). Assuming, half of ferry passengers are 
visitors to Zanzibar^^, the annual number of infectious visitors varied 
from 8.5 to 8.8 per 1,000. The net visitor importation rate is weighed 
by the visitor length of stay respectively to the duration of infection 
and varies between 0.7 and 0.9 infections per 1,000 Zanzibaris 
(Table 1). The total net importation rate to Zanzibar varies between 
1.6 if estimates are based on the risk measured in the main cities of 



each region of origin to 6.7 if estimates are based on transmission risk 
measured in the region (Table 1). These estimates would also vary up 
or down depending on the ratio of residents to visitors. As most 
travel is assumed to occur between Zanzibar and main cities on 
mainland Tanzania, where malaria risk is low, the importation net 
rate of 1.6 infections per 1,000 inhabitants per year is kept as a 
baseline reference for further comparisons. 

In urban Unguja, where an epidemiological survey estimated a 
parasite rate of 0% (95% Confidence Interval: 0-0.2%)'\ the upper 
bound for the net malaria importation rate, if were equal to 0, was 
3.7 per 1,000 inhabitants per year; the lower bound was 0 (Figure 4a). 
The midpoint of the estimates, 1.6 infections per 1,000 inhabitants 
per year, is consistent with the findings from mobfle phone data 
analysis. Assuming malaria importation rates tend to be highest in 
the urban areas that are near the ferry landing, malaria prevalence in 
rural Unguja and Pemba likely provides more information about 
local transmission than about imported malaria. Using the estimate 
from the mobile phone survey of 1.6 imported infections, per 
1,000 per year, R^ ranges between 0 and 0.56 in urban Unguja 
(Figure 4a), between 0.71 and 0.91 in rural Unguja (Figure 4b), 
and between 0.92 and 0.98 in Pemba (Figure 4c). 

Relationships between vulnerability, receptivity and malaria 
prevalence in Zanzibar. The models outlined here iflustrate the 
relationship between the net importation rate, Rc and the pre- 
valence rate for contexts outside of Zanzibar. For the range of 
malaria importation rates (0 to 10 incoming infections per 1,000 
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Figure 3 | Relative importance of residents on malaria importation. 

Contribution of resident travelers to the net malaria importation rate as a 
function of the number of imported infections from visitors versus 
residents (SJSy.) and the length of stay for visitors ( Ty). The solid red line 
represents the contribution threshold below which residents contribute to 
more than 50% to malaria importation. 

per year) and Rc {0 to 1) values used here, parasite rate estimates are 
found to be below 8% (Figure 5a). Assuming the malaria importation 
rate to Zanzibar is less than 1.6 cases, malaria prevalence would not 
exceed 2% unless Rc was close or equal to 1 (Figure 5a). 

Imported malaria not only brings malaria parasites to the island, 
but it also adds to local transmission^^ As Rc and the net importation 
rate increase, the proportion of local infection amid all infections 
increases. Under our current estimation of net importation rate to 
Zanzibar of 1.6 and as Rc is close to one in rural Zanzibar, even low 
undetected net importation rate triggers important local non- 
endemic transmission with high parasite rate that might delay mal- 
aria elimination (Figure 5b). Nevertheless, the overall estimates of Rc 
are below 1, which suggests that if malaria importation ceased. 



endemic malaria transmission would eventually end except in a 
few residual local transmission foci. 



Discussion 

Strategic priorities of national malaria control programs vary 
depending on many factors, including the natural potential for 
endemic malaria transmission, the levels of control that have been 
achieved, and the proximity to other malaria endemic countries. In 
countries where malaria is rare, malaria importation is of primary 
concern. Countries with well-developed health systems, good sur- 
veillance through national case reporting, and no local malaria can 
simply count the number of malaria cases that occur. In the United 
States in 2008, 1,298 detected malaria cases were imported mainly in 
tourists inadequately protected by chemoprophylaxis and immi- 
grants, and none triggered any further transmission^'*; autochthon- 
ous cases do occur, but only rarely^^. Countries that have recently 
reduced their malaria burden, but that have strong human travel 
connections with other malaria endemic countries^^ with poorly 
developed health systems, or some level of ongoing transmission 
must devise alternative methods for estimating malaria importation 
rates. Key information is required to set up appropriate measures: 
contributions of various routes of malaria importation, reliable esti- 
mates of malaria importation rates and of the potential for local 
transmission, and the total impact of these quantities on the local 
endemic and non-endemic transmission. These estimates must also 
be made in a highly dynamic environment as neighboring countries 
also control malaria. Our analysis suggests that citizens of Zanzibar 
who travel to malaria endemic regions are likely to be the single most 
important source of importation. Mobile phone data analysis esti- 
mates that less than 1.6 malaria infections are imported per 1,000 
Zanzibar residents every year, with a credible value of less than 3.7 
infections per 1,000 people per year according to mathematical mod- 
els. This analysis suggests that Rc is below 1 in most places, as most 
local transmission is driven by imported malaria. Therefore, our 
results suggest that, if current levels of control are sustained, malaria 
would eventually disappear from most of the island, except possibly 
in a few residual foci if no imported infections occurred. 

Therefore, enhanced prevention and control strategies should aim 
at mitigating pathogen importation to move towards malaria elim- 
ination and prevent reintroduction. The government could provide 
free chemoprophylaxis through the health system to Zanzibar resi- 
dents travelling to endemic mainland using drugs that minimize the 
risk of developing resistance^^. As immunity to malaria parasite 



Table 1 | Quantitative assessment of imported malaria from mainland Tanzania (Five Zantel regions: Arusha, Dar es Salam, Dodoma, 
Mbeya a nd Mv/a nza ) to Zanzibar (Pemba and Unguja).^rOnd (3^ respectively represents the total number of imported infections brought to 
Zanzibar by residents or visitors infected on mainland Tanzania (importation rate) with T^the median time spent by residents on mainland 
Tanzania and Ty the median time spent by visitors on Zanzibar; Ny represents the total number of visitors travelling monthly to Zanzibar; 
3^ represents visitor net importation rate, the number of infected visitors weighed by the duration of the trip, Ty, relatively to the duration 
of infection. The risk of infection is calculated with risk estimates from the main city or from the rest of the region of origin on mainland 
Tanzania. 
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Figure 4 | Quantitative assessment of malaria importation. Vulnerability (the number of imported malaria per year per 1,000 inhabitants, S) and 
receptivity (the effective reproductive number, R^) in three different epidemiological settings in Zanzibar: (a) urban Unguja with a parasite rate of 0% 
(95% Confidence Interval: 0-0.2%), (b) rural Unguja with a parasite rate of 0.5% (0.3-0.9), and (c) Pemba with a parasite rate of 1.5% (1-2.1). The solid 
black line represents the mean annual number of imported malaria cases and the dashed black lines the upper and lower bounds of the 95% confidence 
interval represented by the grey area. The red dotted line represents the empirical estimation of 3 in Zanzibar according to travel assessment survey, which 
is of 1.6 imported malaria cases per 1,000 inhabitants per year. 



Local 

a) PR(%) u\ infection 

^ (%) 




Figure 5 | Quantitative assessment of the impact of malaria importation on local malaria dynamics. The effect of vulnerability (the number of imported 
malaria per year per 1,000 inhabitants, 3) and of receptivity (the effective reproductive number, R^) on, (a) the parasite rate (PR), and (b) the proportion 
of local infection among all malaria infection. The red dashed line represents the empirical estimation of 3 in Zanzibar according to the travel assessment 
survey, which is of 1.6 imported malaria cases per 1,000 inhabitants per year. The solid red line represents the contribution threshold above which local 
infections contribute to more than 50% to total malaria prevalence 
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declines in low endemic settings, such travelers are even more vul- 
nerable to infection, stressing the rationale to distribute prophylaxis. 
In addition, each case of infection should be reported to health 
authorities, which highlights the need for enhanced surveillance^^. 
Measures should include border screening especially at ports of land- 
ing, active case detection and investigation to set up appropriate 
reactive measures in the environment of the cases. Surveillance tools 
rely on the availability of sensitive and specific diagnostic tools such 
as rapid diagnostic tests (RDTs), which require further improve- 
ment^^. Finally, vector control measures could mitigate the importa- 
tion of infectious anopheles, for instance by using insecticide in 
planes flying from malaria endemic region^°, even though this is 
likely to be a minor issue. Some of these measures may be technically, 
economically and operationally challenging and require prior feas- 
ibility assessment combined with public health education campaigns. 

Imported malaria contributes differently to local dynamics 
depending on the number of incoming pathogens and the length 
of stay as addressed here, but also on when and where malaria is 
imported. Formal and informal human movements are increasingly 
becoming amenable to study as humans use mobile phones. 
Nevertheless, mobile phone usage data are subject to biases based 
on data availability, ownership and usage patterns. Principal malaria 
transmission season on mainland starts between November and 
January and ends approximately between April and June. 
Therefore, the current survey extending over a three-month period 
from October to December only captures a subset of potential mal- 
aria importation. Furthermore, the available Zantel mobile records 
did not provide further information on local cell tower usage, so it 
was not known where residents and visitors stayed in Zanzibar. As a 
result, even though malaria importation could be quantified for the 
main regions of Zanzibar, no local foci of residual transmission can 
be identified at higher spatial resolutions. Finally, even though gen- 
eral population movement to and from Zanzibar mainly rely on ferry 
transportation, informal movement on small fishing boats also 
represent a pathway for malaria importation and can only be par- 
tially captured through mobile phone records. These movements 
create hotspots of malaria importation mainly in coastal communit- 
ies. As health authorities know these hotspots, the Zanzibar Malaria 
Control Programme (ZMCP) can intervene and mitigate malaria 
importation in these communities. 

Some of these limitations have been addressed by cross-validating 
estimates with two completely independent methods using different 
data. Despite the uncertainties, mobile phone usage records have 
provided useful and reasonably accurate information about the 
lengths of stay in Zanzibar and the relative importance of various 
routes of malaria importation. These estimates could be improved 
through further studies focusing on the places where people stay in 
Zanzibar, and gathering mobile phone data over a longer study per- 
iod to better take into account malaria seasonality. To move towards 
elimination, additional tools should be used to complement current 
assessment especially in terms of risk mapping. The Malaria Atlas 
Project (MAP) provides with endemicity maps based on prevalence 
surveys at a global, continental or country level. Using mathematical 
models^^"^^, these maps were converted to EIR risk maps, a more 
relevant measure for risk transmission^^. Nevertheless, as Zanzibar 
is getting closer to elimination, finer scale mapping based on 
environmental, epidemiological, population, vector, social and eco- 
nomical factors is required to identify residual foci in order to imple- 
ment rapid localized response. 

As very low levels of transmission are reached, programs move 
from a focus on control to a focus on pre-elimination and elimina- 
tion, and finally prevention of re-introduction. Quantitative 
approaches to human movement and the spread of malaria are 
increasingly useful for government to appropriately decide when 
shifting efforts to stemming importation. Mobility patterns can be 
used to estimate the import of infectious disease into geographic 



regions with low endemic states. The importance of control methods 
will depend on the vectors and immuno- epidemiology of the patho- 
gen, as for dengue which has a shorter infectious period, daytime 
biting vectors, and life-long immunity^^. As methods improve for 
quantifying malaria importation and transmission, they can be 
applied more widely to advise countries on strategic planning across 
the region and worldwide. 

Methods 

Contributions of residents, visitors and immigrants to imported malaria. The 

standardized contributions of imported malaria from all four routes can be added to 
estimate the malaria importation rate, denoted d, reported as the equivalent number 
of untreated new human infections, per thousand populations, per year. 
Quantitatively, the way residents or visitors influence malaria dynamics mainly 
depends on the number of infections brought by residents and visitors {dr and dy), and 
the length of stay of visitors ( Ty) relative to the duration of infections in residents ( 1/r) . 
Assuming infectivity is constant over the duration of infection, the relative 
contribution of one imported malaria case by a resident compared to one by a visitor 

Sr/r 

is estimated according to the following formula: — — — . 

3r/r-\-SyTy 

Estimating the vulnerability of Zanzibar. Estimating the net malaria importation 
rate from mobile phone usage data. Anonymized phone call record data, which 
included both the time each call was made and the location of the nearest mast that 
each call was routed through, was used to construct trajectories of the movements of 
individuals over time^^. About 9 out 10 Tanzanians have access to a mobile phone^^' 
and the Zanzibar telecom (Zantel) mobile phone operator has approximately a 99% 
market share on Zanzibar^^, this likely comprises a representative sample of Zanzibari 
residents. This does not hold for Tanzania as a whole, as Zantel had only 
approximately a 10% share of the market on the mainland during the data period. 
Mobile phone usage records encompassing a three-month period from October to 
December 2008 were obtained, as the company only archives three months of records. 
This covers the busiest period in terms of travel to and from Zanzibar, and therefore 
enables conservative upper estimates of malaria importation. Zanzibar residents were 
identified based on the assumptions that they made the majority of their calls on 
Zanzibar; the rest were considered to be visitors. Each individual call and message was 
spatially referenced to one of six areas: Arusha, Dar es Salaam, Dodoma, Mbeya, 
Mwanza and Zanzibar (Figure 1). Only those making regular calls (more than four a 
month) with apparent overnights trips were considered, as day trips are of little risk 
considering the mosquito night biting habits^^. The final dataset included 67,236 trips 
made by visitors between Zanzibar and mainland and 73,095 trips made by for 
Zanzibari residents, 88% of Zanzibari users did not travel long distance during the 
study period. 

Mobile phone data inform on individual trips made by Zanzibar residents to 
endemic areas on the mainland, more specifically, where they went and how long they 
stayed (length of stay, T^). The level of malaria endemicity can be measured by the 
daily entomological inoculation rate (EIR or £), which corresponds to the number of 
infectious bites received per person per day. Recent research efforts have provided 
spatial quantification of P. falciparum endemicity enabling estimation of EIR at the 
locations where Zanzibar residents visit^*^' Two types of population-weighted 
estimates for EIR were provided based on the endemicity level for each Zantel region 
either for the entire regions of origin, or their principal cities disaggregated by month 
to take into account transmission seasonality. For residents, the individual probability 
of getting infected by malaria parasites during each trip is described by the following 
formula for the whole trip attack rate that takes into account heterogeneous biting 

= 1 — ( 1 + (xbETr) ^^"^^ h represents the probability that an infectious bite leads 
to infection, and a is an index of heterogeneous risk during a trip in an area with 
average daily exposure £, and a= 4.2 so that approximately 20% of the population gets 
80% of all the risk^°. Combining the information collected with mobile phones to the 
formula of heterogeneous biting allow for estimating the total number of imported 
infections from residents, br, over all iV^ people trips made in during the study period, 

which is br= ^Pr. 

As the mobile phone data do not cover all visitors to Zanzibar, the risk from each 
individual visitor trip provided by mobile phone data was extrapolated to the number 
of visitors travelling each month to Zanzibar provided by ferry traffic data. The 
number of infected visitors traveling to Zanzibar from mainland Tanzania, by, is 

calculated according to the following formula b^ = with Ny the overall number 

of visitors traveling to Zanzibar according to ferry traffic data and X the parasite rate 
for each region of origin or their principal cities, disaggregated by month. To compare 
the respective impact of imported malaria from visitors on local transmission 
to the one from residents, the number of infected visitors was standardized using 

the formula b^ = rTyXi, with Ty the duration of their trip, and r approximately 
0.005 d-\ 

Estimating the net malaria importation rate from a mathematical model. The second 
method, which is based on a model of importation and transmission, estimates the 
malaria importation rate as a function of malaria prevalence and Rc. The model 
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considers only imported malaria cases that are not treated. The dynamics of malaria 
infection at low intensity can be described by a slightly modified Ross-Macdonald 
model"^"' At low prevalence, the model is simplified to X = RcXr(l —X) — rX^^. The 
expression reflects the local transmission of the disease. But the overall malaria 
dynamics should also encompass imported malaria and the expression can be 
improved by explicitly adding per-capita importation rate 5: 
X = RcXr(l -X) + (5(1 -X) - rX. At the equilibrium, X = 0, and the net malaria 

importation rate is estimated by the following formula: d = rX ^ — — — . 

Relationships between vulnerability (S), receptivity (Re) and malaria prevalence in 
Zanzibar. The estimates of receptivity and vulnerability influence malaria prevalence 
and the proportion of local infection amid all infections. Malaria prevalence (X) at 
equilibrium is estimated by solving RcXr(l — X) + (5(1 — X) — rX = 0. In our model, 
the local force of infection is described by R^Xr, therefore the proportion of local 
infection amid the total number of infections in the area of interest is described by the 

following formula ^^^^ ^ . 
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